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,The scope of this paper is to perform a detailed spectroscopic study of the northern O-type supergiant HD 195592. We use a large 
Sample of high quality spectra in order to investigate its multiplicity, and to probe the line profile variability. Our analysis reveals a 
clear spectroscopic binary signature in the profile of the He i A 6678 line, pointing to a probable O + B system. We report on low 
amplitude radial velocity variations in every strong absorption line in the blue spectrum of HD 195592. These variations are ruled 
by two time-scales respectively of 5.063 and about 20 days. The former is firmly established, whilst the latter is poorly constrained. 
We report also on a very significant line profile variability of the H/3 line, with time scales strongly related to those of the radial 
velocities. Our results provide significant evidence that HD 195592 is a binary system, with a period that might be the variability 
time-scale of about 5 days. The second time scale may be the signature of an additional star moving along a wider orbit provided its 
mass is low enough, even though direct evidence for the presence of a third star is still lacking. Alternatively, the second time-scale 
may be the signature of a variability intrinsic to the stellar wind of the primary, potentially related to the stellar rotation. 

fi^ey words: stars: early-type, binaries: spectroscopic, stars: individual: HD 195592 
■PACS: 95.85.Kr, 97.10.Me, 97.20.Pm, 97.80.Fk 
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Introduction 



Most O-type stars are found to be part of open clusters, where 
^ they are believed to be formed. However, some stars are charac- 
terized by large peculiar space velocities and seem to have been 
I ejected from their birth place. Such objects - the so-called run- 
away stars - are important in the sense that they provide the op- 
T— I portunity to investigate quite unusual events, namely processes 
OO at the origin of stellar ejection from open clusters. Mainly, 
two scenarios have been proposed to explain the peculiar kine- 
matics of runaway stars. On the one hand, one may consider 
a supernova explosion in a close binary (or higher multiphc- 
■ity) system, re s ulting in the ejection of the secondary dZwickv , 
Il957t iBlaauwl Il961 ). In the context of this scenario, the re- 
sulting ne utron star may remain bound to the secondary. Ac- 
. !^ cording to lPortegies Zwart about 20 - 40 % of the run- 

aways may indeed be binaries harbouring a compact compan- 
ion. On the other hand, the ejection of the star can be the re- 
6ult of dynamical inte ractions in a dense open cluster (see e.g 



Table 1: Observing rans. The first and second columns give the name of the 
campaign as used in the text. The next columns yield the number of spectra 
obtained, the time elapsed between the first and the last spectrum of the run 
(AT), the natural width of a peak of the power spectrum taken as 1/AT, and 
finally the mean signal-to-noise ratio of each data set. 



Obs. run 



N 



AT 
(d) 



AVnat 
(d-1) 



S/N 
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Oct.2004 


4 


9.2 


0.11 


150 


Sum.2005 


7 


8.0 


0.13 


450 


Sept.2006 


20 


4.2 


0.24 


450 


Oct.2006 


2 


2.0 


0.50 


380 


Aut.2007 


48 


27.0 


0.04 


500 



status. They discovered that seven of them have been ejected 
from their birth place. Three of them were also reported to be 
associated with bow shocks due to their interaction with the in- 



Leonard and Duncanl [l990). In the latter case, a binary fraction terstellar medium ( Noriega-Crespo et al. , 1997 ). Among these 



of about 10 % is expected. Multiplicity studies of runaway stars 
are needed in order to check whether the dynamical interaction 
scenario is sufficient to explain the observed binary fraction, or 
if the supernova explosion scenar io is needed. 



In a recent spectroscopic studv. iMcSwain et al.l ( 120071) inves- 
tigated a series of field OB stars in order to check their runaway 
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runaways, one finds H D 195592 (09 . 5Ia) to which this paper is 
devoted. The study of Mason et al. ( 1998') d id not reveal any 
visual companion, whilst Plaskett and Pearcel ('l93l') suggested 
it may be a spectroscopic binary. [McSwain et al. (2007) re- 
ported also on radial velocity variations but the binary status 
of HD 195592 still needed to be confirmed. Our target star is 
located at a distance of about 1 . 1 kpc, and may be originatin g 
from the open cluster NGC 6913 dSchilbach and RoseilEooil) . 

In this paper, we report on an intensive multiplicity study 
of HD 195592 based on a high quality time series of spectra 
obtained in the blue domain. We first present our analysis, in- 
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Table 2: Journal of the observations of HD 195592, The heliocentric Julian date is given as HJD - 2 450 000, and the date (yyyy/mm/dd) is that of the beginning of 
the night. For each observation, the mean radial velocity expressed in km s" Ms given. The lines used to computed the mean radial velocity are He i -1 447 1 , Si m /U 
4552,4567, Cm 4647,4650, and Hei A. 4713. 



No. 


HJD 


Date 


RV 


No. 


HJD 


Date 


RV 


No. 


HJD 


Date 


RV 


1 


3286.301 


2004/10/07 


-36.0 


28 


3983.478 


2006/09/04 


-30.2 


55 


4410.244 


2007/11/05 


-20.4 


2 


3287.386 


2004/10/08 


-28.2 


29 


3984.310 


2006/09/05 


-29.6 


56 


4410.321 


2007/11/05 


-24.0 


3 


3290.450 


2004/10/11 


-26.3 


30 


3984.390 


2006/09/05 


-31.2 


57 


4411.251 


2007/11/06 


-37.3 


4 


3295.455 


2004/10/16 


-21.0 


31 


3984.465 


2006/09/05 


-33.6 


58 


4411.283 


2007/11/06 


-30.6 


5 


3547.578 


2005/06/25 


-16.1 


32 


4033.258 


2006/10/24 


-21.5 


59 


4411.314 


2007/11/06 


-30.8 


6 


3548.548 


2005/06/26 


-31.5 


33 


4035.255 


2006/10/26 


-33.1 


60 


4412.242 


2007/11/07 


-16.7 


7 


3550.573 


2005/06/28 


-26.7 


34 


4396.268 


2007/10/22 


-34.8 


61 


4412.267 


2007/1 1/07 


-21.4 


8 


3551.533 


2005/06/29 


-33.2 


35 


4396.343 


2007/10/22 


-31.3 


62 


4413.242 


2007/11/08 


-20.8 


9 


3553.548 


2005/07/01 


-26.6 


36 


4397.232 


2007/10/23 


-21.0 


63 


4413.330 


2007/11/08 


-23.5 


10 


3554.598 


2005/07/02 


-21.5 


37 


4397.290 


2007/10/23 


-21.0 


64 


4414.264 


IWlll 1/09 


-28.5 


11 


3555.555 


2005/07/03 


-43.6 


38 


4400.296 


2007/10/26 


-33.9 


65 


4414.368 


2007/11/09 


-25.4 


12 


3980.302 


2006/09/01 


-32.4 


39 


4401.243 


2007/10/27 


-29.8 


66 


4415.246 


2007/11/10 


-37.6 


13 


3980.371 


2006/09/01 


-33.7 


40 


4401.295 


2007/10/27 


-24.2 


67 


4415.333 


2007/11/10 


-33.4 


14 


3980.437 


2006/09/01 


-35.4 


41 


4402.235 


2007/10/28 


-26.0 


68 


4415.363 


2007/11/10 


-39.2 


15 


3980.512 


2006/09/01 


-35.4 


42 


4402.289 


2007/10/28 


-19.8 


69 


4416.228 


2007/11/11 


-33.0 


16 


3981.366 


2006/09/02 


-36.8 


43 


4403.276 


2007/10/29 


-18.5 


70 


4416.301 


2007/1 1/1 1 


-27.9 


17 


3981.442 


2006/09/02 


-33.9 


44 


4403.343 


2007/10/29 


-17.7 


71 


4417.223 


2007/1 1/12 


-36.1 


18 


3982.306 


2006/09/03 


-23.3 


45 


4405.245 


2007/10/31 


-28.7 


72 


4417.248 


2007/11/12 


-37.7 


19 


3982.347 


2006/09/03 


-23.5 


46 


4405.302 


2007/10/31 


-29.9 


73 


4418.223 


2007/11/13 


-29.6 


20 


3982.413 


2006/09/03 


-21.6 


47 


4406.236 


2007/11/01 


-20.8 


74 


4419.225 


2007/11/14 


-31.3 


21 


3982.440 


2006/09/03 


-24.2 


48 


4406.291 


2007/11/01 


-18.9 


75 


4419.250 


2007/1 1/14 


-31.5 


22 


3982.496 


2006/09/03 


-23.7 


49 


4407.244 


2007/11/02 


-20.3 


76 


4421.221 


2007/11/16 


-34.3 


23 


3983.330 


2006/09/04 


-26.4 


50 


4407.324 


2007/1 1/02 


-23.9 


77 


4421.305 


2007/11/16 


-30.5 


24 


3983.356 


2006/09/04 


-28.2 


51 


4408.242 


2007/11/03 


-15.5 


78 


4422.230 


2007/11/17 


-15.2 


25 


3983.376 


2006/09/04 


-28.2 


52 


4408.320 


2007/11/03 


-11.9 


79 


4422.327 


2007/11/17 


-13.4 


26 


3983.421 


2006/09/04 


-27.8 


53 


4409.242 


2007/11/04 


-20.2 


80 


4423.225 


2007/11/18 


-20.5 


27 


3983.459 


2006/09/04 


-29.1 


54 


4409.318 


2007/11/04 


-19.2 


81 


4423.292 


2007/11/18 


-23.6 
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eluding radial velocity measurements and variability tests (Sec- 
tion|3]l. We then discuss the results and present our conclusions 
respectively in Sections|4]and|5] 

2. Observations and data reduction 

Spectroscopic observations were collected at the Observa- 
toire de Haute-Provence (OHP, France) during several observ- 
ing runs from October 2004 to November 2007. All spectra 
were obtai ned with the Aurelie spectrograph fed by the 1 .52 m 
telescope ( Gillet etaD . [l994) . The typical exposure time was 
about 20 - 30 minutes, depending on the weather conditions. 
Our collection of blue spectra (between 4460 and 4890 A ) is 
described in Table[T| All o ur data have been treated foll owing 
the procedure described by iRauw and De Becken (120041) . The 
signal-to-noise ratio of our spectra was estimated in a region 
devoid of spectral lines, between 4790 and 4800 A . We note 
that the quality of the spectra is rather uniform within each ob- 
serving run, with a dispersion in signal-to-noise ratio lower than 
about 20 % with respect to the mean value given in Table[T] 

The complete data set in the blue contains 8 1 spectra, spread 
over a time interval of about 3 years. We managed the obser- 
vational strategy to cover many time-scales, from a few hours 
up to a few years. In 2006, we obtained indeed several spectra 
per night during a few nights to cover short time-scales. In the 
2007 observing run, we obtained generally 2 spectra per night 
during almost 28 consecutive nights. The complete journal of 
observations is given in Table|2] 

In addition, we also obtained 4 spectra in the red domain 
(between 6340 and 6780 A) using the same setup. In this spec- 
tral domain, the resolving power is of the order of 1 1000. The 
HJD - 2450000 at mid-exposure are respectively 3549.571, 
3555.377, 3555.523, and 3652.417. 



3. Results 

The mean blue spectrum of HD 195592, i.e. between 4460 
and 4890 A , with line identification^ is shown in the upper 
panel of Figure[T] We first searched for a significant variabil- 
ity in the spectral lines using the Tem poral Variance Spect rum 
(TVS) analysis method developed by Fullerton et al. ( 1996h . In 
the lower panel of Figure[Tl we show that we detect a signifi- 
cant variability (above the 99 % confidence level) for most of 
the lines, with a higher TVS amplitude for the strongest lines. 
On the other hand, the red spectrum of HD 1 95592 is dominated 
by the emission of H a (left part of Figure|2]l, and presents also 
Nil /1 66 10 in emission and He 1/16678 in absorption. The in- 
spection of our sample of red spectra revealed immediately an 
asymmetric profile for the latter line (right part of Figure|2]i, 
suggesting a partial deblending of two lines reminiscent of the 
case of a double-line spectroscopic binary. 

The temporal analysis of our blue spectral series was per- 
formed using the generalized Fourier method derived by 



'We note the presence of two weak unusual em ission lines respectively at 
4788 and 4803 A . We attribute them to Nii (see e.g. lKa\i[T96^ . 



Heck et all il9S5), and revised by lOosset et al] (l200lh . This 
technique is especially adapted to the case of unequally spaced 
data. This method has been applied to the radial velocity time 
series (see Section lTTT l. and to line profiles (see Section [T2] l. In 
the latter case, our approach consisted in computing the power 
spectrum at each wavelength step leading to a two-dimensioned 
power spectrum. Mean periodograms were then computed to 
obtain an overview of the power as a function of frequency 
for the whole wavelength interval considered in the temporal 
analysis. We also applied prewhitening techniques considering 
candidate frequencies in order to check their capability to re- 
produ ce the spectral variations (see e.g. iDe Becker and Rauw . 
'2004', for applications of these techniques to line profile vari- 
ability analyses). 

3.1. Radial velocity time series 

We measured radial velocities on most of the prominent ab- 
sorption lines present in our blue spectra (see upper panel in 
Figure©, namely He i /I 4471, Sim AA 4552,4567, Cm AA 
4647,4650, and Hei A 4713, by fitting Gaussians. We applied 
the generalized Fourier method to radial velocity time series for 
each individual line, and for the mean of the radial velocities 
obtained for all these lines (see Table|2]i. The power spectra are 
clearly dominated by two peaks, a first one at a frequency (vi) 
of 0.1975 d ' and the second one (V2) close to 0.05 d"' (along 
with their 1 -v aliases). Our time series allowed us to investigate 
frequencies up to about 10 d~^ However, as the strongest peaks 
were located in the frequency interval between and 1 d we 
will restrict the discussion to time scales larger than 1 d. The 
lack of significant power at frequencies larger than Id' is in 
agreement with the absence of detected variations in time series 
made of spectra collected during a same night (see the journal 
of observations in Table|2] for details on the time coverage of 
our data). A summary of the main frequencies reported for the 
radial velocity time series is given in Table|3] These two fre- 
quencies correspond respectively to time-scales of 5.063 and 
about 20 d (respectively Pi and P2 hereafter). We also applied 
the simultaneou s searc h for periods method used for instance 
by iOosset et al.l (l200lh . This method constitutes some kind 
of higher order Fourier technique, taking into account the ex- 
isting correlations between different frequencies. As a result, 
we found again the two frequencies v\ and V2 reported above. 
We note that the latter approach seems to priviledge a value of 
0.055 d ' for V2 even for the mean radial velocity time series, 
although the sequential frequency removal approach pointed to 
slightly different values. 

In order to check for the plausibility of these frequencies, we 
randomly removed about 25 % of the data from our radial veloc- 
ity time series and we repeated the Fourier analysis. We found 
that the main peaks corresponding to v\ and vi are still clearly 
present in the power spectrum. The latter result lends strong 
support to the idea that the presence of these peaks is not related 
to a particular realization of a stochastic process unrelated to an 
actual physical behaviour. As the simultaneous period removal 
approach is more reliable than the sequential one for the study 
of multiperiodic processes, we will favor V2 = 0.0550 d"' (cor- 
responding to P2 = 18.182d) in the following discussion, even 
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Figure 1: Blue spectrum of HD 195592 between 4460 and 4890 A . Upper part: mean spectrum with line identifications. Lower part: TVS with its 99 % confidence 
level illustrated by the horizontal dashed Une. 




Figure 2: Left part: Red spectrum of HD 195592 between 6540 and 6100 A . Right part: Line profile of Hei /i 6678 A respectively at HJD 2 453 459.571 (top) and 
HJD 2 453 652.417 (bottom). The line shape of the lower profile strongly suggests a partial deblending of two spectral lines. 
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Figure 3: Power spectra of the He: A 4471 (upper panel) and mean (lower 
panel) radial velocity time series. Both power spectra are dominated by vi and 
V2, along with their 1 - v aliases. The frequency is expressed in d"' and the 
amplitude in km s"' . 



Table 3: Main frequencies expressed in d"' derived from the Fourier analysis 
of the radial velocity time series. The frequencies marked with a * correspond 
to peaks slightly weaker than that of their 1 - v alias. 



Spectral line 


Vl 


V2 








He I A 4471 


0.1975* 


0.0415 


Si III A 4552 


0.1975 


0.0485 


Si III A 4567 


0.1975 


0.0440 


C III A 4647 


0.1975 


0.0475 


C III A 4650 


0.1975 


0.0475 


He I A 4713 


0.1975 


0.0415 


Mean 


0.1975 


0.0440* 



though this value deviates significantly from those quoted in 
Table[3] 

We fitted sine functions with the period Pi to our radial ve- 
locity time series. We then subtracted the fitted radial velocities 
from the measured ones, and we folded the residual radial ve- 
locities with the period P2. We then proceeded the same way 
but by fitting the measured radial velocities with a period P2 
sine function, and folding the residual radial velocities with Pi . 
We repeated the same procedure for the measurements obtained 
from every line, and for the mean of all lines used for the mea- 
surements. The resulting residual radial velocity curves sug- 
gest trends a priori not compatible with noise. Even though 
the semi-amplitude of the radial velocity variations are of the 
same order of magnitude than the expected typical error on 
our measurement^ these curves suggest that two simultane- 



Figure 4: Fourier analysis of the H /3 line betwen 4860 and 4863 A . Upper 
panel: mean power spectrum between and 5 d"' . Middle panel: mean power 
spectrum obtained after prewhitening using the frequency of the highest peaks 
of the initial power spectrum, i.e. 0.1975 and 0.15375 d"' . Lower panel: spec- 
tral window directly related to the sampling of the time series. 



ous periodic phenomena are responsible for the variations re- 
vealed by the temporal analysis. The consistency of the results 
obtained for the different lines concerning v\ indicates that its 
presence is firmly established by our temporal analysis. How- 
ever, the deviations observed from one line to the other for V2 
deserve some caution concerning the corresponding variability 
time scale. The different values close to V2 reported in Table[3] 
allow a time scale ranging between about 20 and 24 days, down 
to 18 days according to the results of the simultaneous search 
for periods. 



3.2. Line profile variability 

Beside the radial velocity variations discussed above, we also 
detect a substantial line profile variability in H j3. Our tempora l 
analysis using the same tools as De Becker and Rauwl (120041) 
reveals a variability dominated by a frequency of 0. 15625 d~' 
(v3, corresponding to a time scale close to 6.4 d, hereafter 
P3). After prewhitening with this frequency, the residual power 
spectrum is dominated by vi. When the power spectrum is 
prewhitened using v\, V3 is slightly shifted to a lower value, i.e. 
0. 15375 d ' . The resulting periodogram after subtraction of the 
two reported time scales is presented in Fig. ID Considering the 
uncertainties on the frequencies emerging from our temporal 
analysis, V3 is compatible with the difference V2 - v\ (mostly 
if V2 is close to 0.044 d"'. In that case, the variations in the 
line profile of H/? would therefore be totally determined by the 
frequencies v\ and V2. However, the prewhitening using v\ and 
V2 gives a residual power spectrum dominated by a peak close 



-We estimate the peak-to-peak amplitude for the variations of the radial 



velocities of Diffuse Interstellar Bands (DIBs) to be of the order of 6 km s 



to V3, with a varying amplitude depending on the chosen value 
for V2. Alternatively, V3 may thus be related to an independent 
process. 

The Fourier analysis of the He i. Si iii and C iii line profiles 
reveals also the presence of strong peaks at frequencies at - or 
very close to - vi and V2- We report also on the presence of 
these frequencies in the central part of He 11 A 4686 (displaying 
a weak central emission surrounded by two shallow absorption 
parts), with a weak residual peak at a frequency close to V1+V2. 
For the red absorption part of the same line, the situation is 
more confused but vi is present. 

H/3 is the only hne in the blue spectrum of HD 195592 that 
presents a strong line profile variability. To illustrate it, we 
have plotted a sample of the profiles in Figure|5] In order to 
check for any trend related to vi, i.e. the most clearly present 
frequency arising from our temporal analysis, we have folded 
the profiles using the corresponding period Pi . The selected 
spectra have been chosen in order to sample as uniformly as 
possible the complete phase interval between and 1. Even 
though some consecutive profiles display obvious similarities 
in the plotted sequence, therefore suggesting a progressive and 
continuous profile variation, it should also be noted that anoma- 
lies are observed as compared to an expected smooth variation 
ruled by a unique and well determined time scale. This may be 
due to the simultaneous impact of P2 on the line profile varia- 
tions. 

As in the case of the radial velocity time series, vi appears 
much more firmly established than V2, as it is found invariably 
in all lines investigated in the blue spectrum, with no deviation 
with respect to the value given in Table[3] 

4. Discussion 

The main result reported above is the signature of a spec- 
troscpic binary with at least two simultaneous variability time- 
scales for the radial velocities. Considering the apparent lumi- 
nosity ratio suggested by the two blended components in the 
right part of Figure|2] the secondary might be a B-type star. As 
the SB2 signature is found in only one line among those acces- 
sible in our data, it is not possible to speculate on the accurate 
spectral classification of the secondary. Both spectral types and 
luminosity classes affect the luminosity ratios of the two stars of 
the system, and it is not possible to translate an equivalent width 
ratio of one single line of the two components into a spectral 
classification for both stars. The fact that asymmetric profiles 
are not obvious in the blue domain points to low amplitudes ra- 
dial velocity excursions (in agreement with the measurements 
discussed in Section lTTl ). In the context of a binary scenario, 
one of the two time-scales (Pi or P2) should be the orbital pe- 
riod. The fact that Pi is much more firmly established than P2 
leads us to consider the former as the better candidate for the bi- 
nary period. However, the origin of the second time-scale needs 
to be investigated. Potential scenarios are considered below. 
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Figure 5 : Plot of a sample of 20 HyS profiles (between 4845 and 4870 A ) folded 
in phase according to the time scale Pi . The phase is specified in each case. The 
To was arbitrarily chosen as the time of the first spectrum of the time series. For 
the sake of clarity, the plots are shifted in normalized flux by an arbitrary value 
of 0.15 units. The spectra have been selected in order to reach a quasi uniform 
phase coverage of the candidate period Pi. The apparently constant absorption 
feature at about 4867A is attributed to N in. 



4.1. The triple system scenario 

A biperiodicity in the radial velocity variations is reminiscent 
of what one may expect from a triple system. In the context of 



6 



this scenario, we can envisage two separate situations: 

Case 1: only the primary spectrum is detected because the 
lines of the other two stars are too weak, or even absent. 
In this case, the observed radial velocity time series is the 
result of two superimposed orbital motions. 
Case 2: the observed line profiles are the result of the blend 
of two or three components associated with the stars of the 
multiple system. 

4.1.1. Case 1 

As a first approximation, we may consider that the observed 
radial velocity variations are due to one of the two orbits, per- 
turbed by the second one. The measured radial velocities are 
therefore representative of the motion of the primary star, and 
Kepler's laws should apply. 

We followed an iterative orbital disentangling procedure in 
order to separate the orbital elements of the two putative orbits. 
Let us define the following quantities, where / goes from 1 to 

aS 

Di : radial velocity measurements 
Ri j : theoretical radial velocities of Orbit 1 at iteration j 
Qi j : theoretical radial velocities of Orbit 2 at iteration j 
Xi j = Dj - Qi j : radial velocities corrected for Orbit 2 at 
iteration j, i.e. approached radial velocities of Orbit 1 
Yij = Di - Ri j : radial velocities corrected for Orbit 1 at 
iteration j, i.e. approached radial velocities of Orbit 2 

As a starting point of the procedure, we estimated the /?, o on 
the basis of a simple sine function with a period equal to Pi. 
We then calculated 7, i, and we computed the first approached 
orbital elements of Orbit 2, usin g the Liege Orbital Solution 
Pa ckage (LOSP, seelSana & Goss et (2009)), as already apphed 
by iDe Becker et al. I (l2006h . On the basis of the latter orbital 
solution, we derived Qn and the corresponding We used 
LOSP iteratively and convergence was generally obtained for j 
= 2 or 3. 

The semi-amplitude of the radial velocity curves of the two 
putative orbits computed by LOSP are respectively equal to 
5.1 +0.8 and 5.5 + 0.8 km s"'. In addition, the two orbits are 
compatible with null or small eccentricities depending on the 
value adopted for P2. The low amplitude of the radial velocity 
variations strongly limits the quality of the results, and the or- 
bital elements are computed from radial velocity series which 
are corrected on the basis of approximated theoretical radial 
velocities carrying inevitably additional uncertainties. The pa- 
rameters we derived request therefore some caution. Finally, 
we note that we obtained a very low value for the projected 
semi-major axis : of the order of 0.5 and 2.0 R© respectively 
for the two orbits. If real, these values either indicate that the 
inclinations of the orbits are very low (i.e. the system is viewed 
pole-on), or that the secondary and tertiary components have 
very low masses compared to the primary. 



4.7.2. Case 2 

When line profiles are made of two blended components re- 
spectively related to two orbiting stars, a slight variation in the 
centroid of a spectral line is expected. In such a situation, the 
measured radial velocities are therefore not a measurement of 
the orbital motion of one star, but would rather represent the 
variations of the centroid of a spectral line made of at least 
two slightly shifted components. In this context, as the com- 
ponents are not disentangled, the measured radial velocity vari- 
ations are not obeying Kepler's laws, and the orbital elements 
derived above in Case 1 do not make sense anymore. Consider- 
ing the significant detection of a secondary component at least 
in He I /I 6678, this scenario seems indeed plausible. 

We had a very close look at line profiles (mainly He i A 447 1 , 
the strongest absorption line in the blue spectrum) in order to 
seek for asymmetries, and slight anticorrelated width and depth 
variations. We detect marginal variations of the depth of the 
line, and slight asymmetries with the red wing of the absorption 
profile steeper than the blue one. However, such an asymmetric 
profile may result from a weak emission component in the red 
wing of the He i line related to the wind material, and is there- 
fore not necessarily due to a blue-shifted secondary absorption 
component. 

If the second time-scale P2 is indeed related to the presence 
of a third - still undetected - star orbiting along a wider orbit, 
the measured radial velocity variations could be explained by 
the modulation of the blended profile of the central binary (dis- 
cussed a few lines above) by the motion of the putative tertiary 
star. In the context of such a scenario, we would be dealing 
with a slightly hierarchical triple system with a rather low pe- 
riod ratio P2/P1 (~ 4). Considering a m ass of about 30 for 
the 09.51 primary (IMartins et al. , 2005i) . and assuming a mass 



of about 20 Mq for the secondary, we applied th e stability crite- 
rion established by^ Leonard and Dunca 311990') for hierarchical 



triple systemfl Our result is that such a system may be dynam- 
ically stable provided the mass of the third star is not larger than 
about 7M0. This upper limit on the third mass translates into 
a spectral type not hotter than a late-B if it is a main-sequence 
star, or even late-A if it is a supergiant. Consequently, even 
though the period ratio of such a hypothetical hierachical sys- 
tem is quite small, its existence should not completely be ruled 
out provided the mass of the third star is low enough. 

4.2. Stellar pulsations 

Line profile deformations responsible for a slight shift of 
the centroid of a line profile can also be caused by pulsations 
of a single star. However, such phenomena ar e expected to 



occur on time scales of a few hours (see e.g . iKambe et al 



1997: Rauw et al.. 2008: De Becker etal., 2008), much shorter 



than the time scales reported above in the case of HD 195592. 
Non-radial pulsations may give rise to frequencies of a fraction 



^^In the present case, N is equal to 81, the total number of spectra in our 
sample. 



^As our data did not allow to determine the spectral type of the secondary, 
our assumption relied on the fact that the secondary is significantly less lumi- 
nous than the primary, but still detected. We note however that the value of the 
inner mass ratio does not have a strong impact on the present discussion. 
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of a day only if beating or resonance is occurring jKaufer et al 
I2OO6),' but this is a priori unlikely to happen without the 
presence of parent peaks at much higher frequencies. We insist 
on the fact that the power decreases rapidly as a function of 
frequency in our power spectra, lending strong support to the 
idea that the variability time scales are of several days. We 
therefore do not consider non-radial pulsations as a plausible 
explanation for the variations revealed by our spectral time 
series. 



4.3. Stellar rotation 

On the other hand, one may also consider one of the re- 
ported variability time-scales to be the rotation perio d. We 



used the relations (1 to 3) from ' De Becker and Rauwl (120041) 
in order to estimate the expected minimum (constrained by 
the critical rotation velocity) and maximum (constrained by 
the projected rotational velocity) rotation perio ds. To do s o. 



we con sidered stellar parameters as prescribed by lMartins et al 
(l2005h for an 09.5 supergiant (M. = 30Mq, R. =23Ro and 
L, =3.1 X IO^Lq), along with the projected rotational veloc- 
ity determined by Conti and Ebbets ( 1977), namely 60km s"'. 
These values lead to minimum and maximum rotation periods 
respectively of 2.7 and 19.4 d. P], P3 and to some extent P2 are 
therefore compatible with the allowed interval. Therefore, we 
cannot reject a scenario where one of the time-scales might be 
related to the stellar rotation. For instance, the presence of P3 in 
H y6 may point to a rotational modulation of the profile produced 
partly in the wind, with V2 resulting from a beating between vi 
and V3. However, if V3 is indeed a true frequency, we would 
have to explain why it is absent from other lines although vi 
is ubiquitous in absorption lines. Alternatively, P3 in H/? may 
still be related to the stellar rotation of the primary whatever the 
origin of P2. 

4.4. HD 195592 as a runaway 



HD 195592 has been classified by iMcSwain et al.l (120071) 
as a runaway on the basis of a peculiar space velocity (Vpec) 
criterion, i.e. Ivp^d > 30 + a-ypec where a-yp^c is the er- 
ror on Vpec- They determined a peculiar space velocity of - 
42.3 + 8.2 km s"'. This value was calculated using radial ve- 
locities similar to ours. The runaway classification proposed 



by iMcSwain et al.l (l2007h is therefore not contradicted by our 
measurements. The parent cluste r of HD 1955 92 remains un- 
known, even though the study of Ide Wit et al. I (I2OO4) revealed 
the presence of a nearby cluster 

As mentioned in Section[T] several runaways have already 
been identifi ed as binaries harbour ing a compact companion. 
For instance. iMcSwain et al. (2007) summarize the orbital ele- 
ments of five spectroscopic binaries that belong to the category 
of runaways. In the case of HD 195592, we report above on the 
detection of a secondary that may be a B-type star Even though 
its spectral type is not determined, the presence of its spectral 
signature rules out a degenerate secondary star We may envis- 
age two situations: 



- HD 195592 is an SB2 system, with the second variabil- 
ity time-scale being related to the stellar rotation. Such a 
binary system may have been ejected through dynamical 
interactions in its parent open cluster. 

- HD 195592 is a triple system, made of a close massive bi- 
nary plus a much lower mass unidentified third star. In 
such a situation, we would be dealing with the only run- 
away triple system identified so far. However, our anal- 
ysis did not provide clear evidence for the presence of a 
third star. The triple system scenario should therefore only 
be considered as a working hypothesis for future investi- 
gations devoted to HD 195592. Considering the low am- 
plitude variations of the radial velocities in this system, a 
much higher spectral resolution is needed in order to (i) 
confirm the presence of a third star and (ii) identify its na- 
turll Such an identification is requested if one wants to 
better understand the nature of HD 195592 in the frame- 
work of the study of runaway stars. 

5. Summary and conclusions 

We found strong evidence for the signature of a spectro- 
scopic binary in the He 1/16678 line of the runaway supergiant 
HD 195592. Using unprecedented spectral time series, we in- 
vestigated the radial velocities and the line profile variability 
in the blue part of the optical spectrum. We report on a low 
amplitude variability of the radial velocities with frequencies 
vi and V2 respectively of 0.1975 d"' (Pi = 5.063 d), and close 
to 0.05 d ' (P2~20d), in all prominent absorption lines. We 
also performed a line profile variability analysis of the H y6 line 
which displays complex variations, dominated by a frequency 
close to 0. 15 d"', in addition to vi that is also clearly present. 

We attribute the period Pi to the orbit of a close binary system 
made of an 09.7I primary with a probable B companion. In this 
scenario, the line profile variations of H/? may be influenced 
by the combined effect of the binary motion and of an intrin- 
sic variability of the primary. The second time-scale (P2) may 
be the signature of the presence of a putative third star moving 
along a somewhat wider orbit. Dynamical stability considera- 
tions suggest a triple system with such a low period ratio may 
exist provided the mass of the third star is low enough. Strong 
evidence for the presence of a third star is however still lacking. 
Alternatively, the second time-scale may be related to variations 
intrinsic to the stellar wind of the supergiant primary, possibly 
due to the stellar rotation. We note that in the context of the 
investigation of runaway stars, the case of HD 195592 deserves 
a particular attention in the sense that it is now confirmed as a 



^At this stage, it it interesting to note that our data cannot reject a triple sce- 
nario including a neutron star as a third object. Such a situation may occur in 
the context of the supernova ejection processes, even though it is not clear so 
far whether such a multiple system could survive to the supernova explosion. In 
addition, the presence of the neutron star on the wider orbit of the hierarchical 
triple system would deserve some attention, as a mass segregation is generally 
observed in high multiplicity system. It is so far not sure that dynamical in- 
teractions consecutive to the supernova explosion could lead to some orbital 
exchange within the system. 
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SB2 system, in contrast with other runaway binaries identified 
as SBl systems such as HD 14633 and HD 15137. 

The results of the present study allow us to draw guide- 
lines for future observations. In order to clarify the nature of 
HD 195592, the signature of the SB2 system should be searched 
for in the complete visible spectrum, investigating as many lines 
as possible, and the putative signature of the third star should 
be scrutinized. In addition, observations in the X-rays may pro- 
vide valuable information as well. HD 195592 is included in 
the field of two ROSAT observations (Seq. IDs RS930941N00 
and RS930942N00) but does not seem to be a bright soft X- 
ray emitter. However, more sensitive observations covering a 
wider energy doamin (e.g. with XMM-Newton) are expected to 
be helpful. 
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